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Soft mode near the critical end point in the phase diagram of two-flavor Nambu-Jona-Lasinio 
(NJL) model is investigated within the leading 1/Nc approximation with Nc being the number of 
the colors. It is explicitly shown by studying the spectral function of the scalar channel that the 
relevant soft mode is the scalar density fluctuation, which is coupled with the quark number density, 
while the sigma meson mode stays massive. 
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I. INTRODUCTION 

Existence of critical end point (CEP) in the QCD phase diagram is recently suggested in several works [1-5], and 
its implications to the phenomena in the heavy ion physics are intensively discussed in order to identify its location 
experimentally [6-14]. 

The critical end point is the end point of the phase boundary determined by the first order phase transition. The 
nature of this point has been discussed so far mainly from the viewpoint of the chiral symmetry. The order parameter 
is the scalar quark condensate, {qq) ^ a, which strictly vanishes in the Wigner phase in the chiral limit. At the 
CEP in question the phase transition necessarily becomes of second order even with the non-zero current quark mass 
m 7^ 0, and the curvature of the free energy w.r.t. the scalar field should vanish. This fact immediately implies the 
emergence of a gapless excitation in the scalar channel. The existence of this second-order transition, however, is 
accidental in the sense that no symmetry change of the ground state is accompanied with this transition. In fact the 
finite current quark mass, which breaks the chiral symmetry explicitly, results in the non-vanishing scalar condensate 
at any point in the plane of the temperature and the quark chemical potential (T — /i) , in contrast to the case with 
m = 0. Only its fluctuation around the equilibrium value becomes unstable at the CEP. This "accidental" instability 
is not related to any symmetry of the system. 

Even in the chiral symmetric world with m — 0, presence of the tricritical point (TCP), which is the counterpart of 
the CEP, is not necessarily required by the spontaneous breaking of the chiral symmetry. The change of the symmetry 
realization of the system constitutes a boundary line in the phase diagram, on which the first-order or the second-order 
transition occurs. The soft mode for the chiral breaking is known as sigma meson, and it becomes gapless on the 
second-order line so as to recover the chiral symmetry together with the Goldstone pion. The TCP exists as a point 
separating the boundary line into the segments of the first-order and the second-order transitions, depending on the 
parameters T and fi. What kind of singular behavior will appear additionally at this TCP, where the sigma mode is 
already gapless? 

The most natural expectation may be that the density fluctuation becomes soft near the CEP because it is the 
end point of the first order transition which accompanies the density gap caused by the strong attraction between 
the quarks. The aim of this paper is to prove this scenario explicitly within the NJL model [15,16] in the leading 
1/Nc approximation, and to show that the nature of the CEP is rather similar to the liquid-gas critical point than 
the chiral transition. The density fluctuation discussed here will be the collective motion of the particle-hole-like 
excitations. One should note that this density fluctuation has a space-like dispersion relation from the kinematical 
reason. Thus, we will be forced to investigate the mode spectrum with this kinematics. If this is the case, some of the 
implications discussed before in the literature [6,10,13] must be reconsidered in relation to the heavy ion experiments; 
at the CEP the sigma meson, which is the chiral partner of the pion, can remain massive. One should study the 
observables distinctive of the space-like density fluctuations, rather than the particle production processes with the 
time-like kinematics. 

The nature of the phase transition is insensitive to the microscopic structure of the system. The most important 
ingredients which determine the nature of the transition are the dimension of the system, the symmetry of the 
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interactions or the dimension of the order parameter(s), and the range of the interaetions. If the study is extended 
to include the dynamical nature, these universarity classes for the static properties split into subclasses depending on 
the existence of the conservation law(s) and mode-mode coupling(s), and so on [17]. It may be, therefore, perilous to 
utilize (e.g.) the Ising model results in order to speculate the dynamical nature of the QCD phase transition. With 
this caution in mind we employ a concrete model of the quarks in this paper to study the dynamical properties of the 
phase transition at the CEP. 

We use the NJL model here as a simple model which can describe the CEP as well as the chiral phase transition 
of the quark system. One of the advantages of the use of the NJL model is that the bosonic collective modes of the 
system are easily and explicitly calculated. Our main conclusions are expected to hold irrespective of the details of 
the model employed here. We should still keep cautious, however, because this model lacks the long-range nature of 
the color interactions, and the gluonic degrees of freedom at all. In this paper we restrict our discussion only for the 
second order phase transition at the CEP (to 7^ 0). More complete analysis treating the TCP, and the vector and 
entropy density fluctuations, will be reported in the forthcoming paper [18] in the NJL model as well as in the linear 
sigma model. 

This paper is organized as follows. In §11 the behavior of the susceptibilities in the NJL phase transition is discussed. 
In §111 the spectral function of the scalar mode is investigated in some details and the relevant slow mode is explicitly 
identified as the scalar density fluctuation. §IV is devoted to summary. 



II. SUSCEPTIBILITIES IN THE NJL PHASE TRANSITION 



We use the simplest version of the NJL model with two flavors {Nf = 2) [19] whose laglangian is, with obvious 
notations, 

jC = q{i^- m)q + g[{qqf + {qil^r" qf]. (1) 

After introducing the auxiliary fields a and tt for qq and qij^r'^q in a standard way, the system pressure at finite 
temperature T and finite quark chemical potential /x is written, within the mean-field approximation with cr=const 
and 7r=0, as 

1 t rl^k 11 

-P{T,„,m;a) = J - Tln(l - n_) - Tln(l - n+)] - (2) 

where n± = (ef(^^i^') + E = VM^ + k^, M = m — a, and u = 2NfNc is the number of the quark species. 

This expression corresponds to the leading 1/Nc approximation with gNc being 0(1), though we set Nc = 3 in the 
following. From the pressure P{T, fj,, m; a) we can calculate all the thermodynamical quantities. 

The ground state at (T, /j.) satisfies the stationary condition for a, which yields the gap (saddle point) equation 



/ 



d^k M ,^ ^ 1 1 

(l-n_ -n+)-|--— cr = 0. (3) 



(27r)3 ^' 2gy 

We define the curvature of the potential w.r.t. the scalar field at the saddle point 

1 fifkl,, . f (fk M2 /I 



fifkl,^ , f (fk {l-n--n+ , .\ 
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= ^^-x'^l>G (4) 

with = dn±/ dE. The last inequality must be satisfied according to thermodynamical stability. The occurrence of 
the second order phase transition is signaled by the vanishing curvature, K = Q. Furthermore, we notice from Eq. (4) 
together with the gap equation, that K ^ ^ at the second order chiral transition approached from the lower 
temperature or density. As we shall see later, K is proportional to the inverse propagator of the static scalar mode 
in the long wave length limit. 

The susceptibilities Xab = d^P/dadb {a,b = T,^,ni) are very important to characterize the nature of the phase 
transition [20,21]. In our approximation these susceptibilities have a simple common structure as 

y , = y(0)ly(°) +Y^1^ (5) 
X.ab — Aam A.mb ' ^ab 1 \^ I 
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FIG. 1. Calculated phase diagram of the NJL model. The transition line of the second order is drawn in the dashed line. 
The phase boundary of the first order transition constitutes a surface shown by hatch. The cross indicates the tricritical point, 
where three lines of second order transition meet. 



where xi^^ are the non-singular parts of the susceptibilities and equal to those of the free quark gas of mass M. From 
this expression we find that all divergences at the second order phase transition originate from the flatness of the 
potential curvature K —^ 0. One should note here that the mixing of the susceptibilities with the scalar mode xim 
(a = T, n) is proportional to M. Thus, at the chiral phase transition where m = 0, the would-be singular parts of 
Xab {o.,b = T,ii) involving the factor l/K, arc less singular than Xrnm- Actually they remain finite (vanish) if the 
transition point {Td^c) is approached from the broken (symmetric) phase because the constituent quark mass M 
vanishes on the boundary and in the symmetric phase. On the other hand, at the CEP (m 7^ 0), M stays finite and 
the critical behavior of these susceptibilities is essentially the same as Xmm- This argument is consistent with the 
analysis based on the general properties of the Landau free energy with the tricritical point [12], as it should be. 

The discussion can be extended a little bit further by changing the thermodynamical variables. Let us choose the 
quark number density pq as an independent variable instead of p,. Then the scalar susceptibility with fixed quark 
number density is related to the susceptibilities given in Eq. (5) via 

Ammy-'- 1 Pq) — yXmmXfJ,lJ. Xm/j.) ~ (g) (o) i (0) " ^ ' 



We find that the K ^ term is canceled out in the numerator. At the chiral phase transition we already pointed out 
that Xfifi ~^ constant, and then the scalar susceptibility with the fixed quark number density Xmm{T, Pq) diverges in 
the same way as Xmm- There the vector susceptibility is decoupled from the scalar one as Xm)* ~ M ^ 0. In case 

of the CEP which we are concerned about, however, Xm^ remains finite, and the denominator and the numerator in 
Eq. (6) diverges in the same way as l/K, which results in the finite scalar susceptibility. This demonstrates that the 
coupling with the number density fluctuation is essential for the divergence accompanied by the CEP. Similarly the 
thermal susceptibility with fixed quark number density, XTT{T,Pq) = (xttXmm ~ X^t)/Xmm' '^^^^ diverge at the 
CEP, but xtt in Eq. (5) does. 

This situation is completely analogous to the thermal susceptibility or the specific heat in the liquid -gas phase 
transition, where the density gap between two phases is identified as the order parameter. The specific heat at 
constant pressure and the isothermal compressibility diverge like \T — Tc\~^ at the critical point. Once we keep the 
number (entropy) density fixed, however, the specific heat at constant volume (the isentropic compressibility) remains 
finite or diverges only weakly like \T — Tc|~") [22]. 

The NJL model has two fundamental parameters: the coupling constant g and the cutoff A. The "physical values" 
of these parameters are usually fixed so as to reproduce the chiral quark condensate {qq) and the pion decay constant 
/tt- Taking the three-momentum cutoff scheme, the critical value gcA^ for the spontaneous symmetry breaking in the 
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vacuum is readily found as cjcA^ = ^f- = 1.64- • • in the chiral limit (m = 0). There is also a critical value for the 
appearance of the TCP, which value is obtained as qtcp^^ = i-e-a = 1-90 • • • from the condition that the err- 
and cr^ -terms simultaneously vanish in the pressure -P(T, /i, 0; a). It is interesting to note that the TCP must appear 
with the finite chemical potential (T/A, fi/A) = (0, 1/e) at g = grcp and that it approaches the T-axis asymptotically 
as g getting larger; this fact proves that the TCP is only possible at finite number density in our model. 

In this paper, for demonstration, we set gA^ = 2.5, and all other dimensionful quantities will be measured in the 
unit of A. The phase diagram of the NJL model in the mean field approximation is shown in Fig. 1. This shares the 
well-known structure of the phase diagram possessing a tricritical point [23] . The character of this diagram around the 
TCP is studied using the Landau free energy and the Schwinger-Dyson equation of the ladder QCD in Ref. [12]. Their 
results on the quark mass dependence of the boundary curves and the behavior of the susceptibilities are consistent 
with those in NJL calculations [18], because the symmetry and order parameters of the system are the same. 



III. RESPONSE FUNCTIONS AND MESON SPECTRUM 



In order to clarify the structure of the relevant modes for the phase transition, we should investigate the spectral 
functions [24-27] . Using the linear response theory, we first calculate the response functions in the relevant channels, 
whose imaginary parts give rise to the mode spectra. 

Let us assume artificially small deviations of the chemical potential and the current quark mass, p,, rh, from their 
equilibrium and physical values, respectively. The response functions for these disturbances are calculated as 

Xa6(«a4,q) = n„b(ig4,q) +n„„(ig4,q)- -2glimb{iqi, fl) , {a,b = n,m). (7) 

Here the polarizations are defined with the imaginary-time quark propagator S as 

j^TY,tricDS{k)S{k-q), (8) 

/d'^k v-^ 
j^-^T22trfcDS{k)ijiS{k-q)iji, (9) 

/(Pk 
j^T J2 tri,uS{k)S{k - q)i^i, (10) 

where = 2lirT (Z e Z ) is the Matsubara frequancy of the bosonic mode, k = (k, k4+iiJ,) = (k, —ujn+ilJ') is the quark 
loop momentum with tt;„ = (2n + l)7rT (n e Z ), and the traces are taken over the flavor, color and Dirac indices. 
The real-time response function is obtained by an analytic continuation, which is simply achieved by a replacement 
iq4 qo+ie after the Matsubara sum. These response functions in the static case should reduce to the susceptibilities 
in Eq. (5). From this fact we find that the static polarization 1106(0, q ^ 0) are the free susceptibilities xfbi^ ~^ 0)> 
which can be also confirmed by explicit calculations. 

The common denominator, ^ — nmm(9ojQ)i is proportional to the inverse of the scalar response function, and 
its static value coincides with the potential curvature K in the long wave length limit. Diagrammatically, this 
denominator corresponds to the sum of the bubble diagrams in the scalar channel. The gapless excitation mode which 
is responsible for the appearance of the CEP and for all divergences of the scalar, vector and entropy fluctuations 
through the coupling, must be contained in the spectrum of this bubble diagram. 

The only task we have to do, therefore, is to investigate the spectrum of the scalar excitation determined by the 
condition 

X^nL « ^ - ^mmiqo, <i) = ^ ' Ji<l) + (4^' ' Q^)I{Qo, q) = 0, (11) 



where 



cPk ^ 1 



f d:'k 



(k^ + M^){{k - qf + A^P) 

1-1 — T7ri rj.-i — T7,r. 1 — n t — rj.r. 

(12) 
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FIG. 2. Sigma (x) and pion (*) masses Ma,Tv at the CEP as functions of the current quark mass m. The constituent quark 
mass M (+) is displayed as well. 



and 



- n_i 



2Ei 



_n+i ^ l_ 



n-2 - n+2 



2Eo 



27r2|q| 



dkk 



+ Tln 



(1 



(1 + e-('^">'°-'")/'^)(l + e-('^">'»+'')/'r) 



(13) 



with £i,2 = v'A^^ + (k±q/2)2, n±i,2 = n±(£;i,2) and 



= ^ M2 + + q2/4 ± |q| |k| . In these equations we 



shifted the three momentum k to make the integral more symmetric and chose a simple convention of the cutoff for 

the mode with finite momentum q as fc^^^ = — q^/4. The stable mode should be obtained as a real solution of 
Eq. (11). The unstable, resonance mode will be found as a complex-valued solution in the unphysical lower half plane 
of go- 

The function /(go, q) is elementary in the thermal field theory [28]. By noting 1 — n+ — n_ = (1 — n_)(l — n+) — 
(pair creation — annihilation) and n_i — ri_2 = (1 — n_2)n_i — (1 — n_i)n_2 (emission — absorption) and so on, we 
find that the first and last terms in Eq. (12) are contributions of the time-like spectrum while the second and third 
come from the space-like spectrum in the spectral representation for /(go, q)- It niay be instructive to note here that 
two limits of the function /(go, q) are different: 



/(0,q^0) 



d^k 1 fl 
(27r)3 4£2 



n+ - 
E 



+ n_ + n' 



(14) 



f (Pk 1 1 - n I - n_ , , 

%o^0,0) = .y^ — 1 . (15) 

The difference is due to the space-like mode contribution; in the first limiting procedure / has the space-like mode 
contribution while in the second case we don't see it. This is the origin of the fact that the massive sigma meson and 
the zero curvature of the free energy can cope with each other [See also Eq. (4)]. 



A. Sigma and pion mcisses 

First we study the excitation energy of the massive mode, which is continuously connected to the sigma meson in 
the vacuum as the parameters (T, ji) are varied. To this end we set g = (go, q) = {lo, 0), following the usual calculation 
of the sigma meson mass in the vacuum. Then J = v J ^^J^.^ (1 — n_ — n+)/E, and together with the gap equation 
we recover the well-known condition to determine the sigma energy, 

__ + (4M2-a;2)/(a;,0) = 0. (16) 

Similarly for the pions, 
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FIG. 3. Calculated spectral function in the scalar mode with (r/A,m/A, |q|/A) = (0.14985,0.01,0.1). (a) Profile at the 
CEP (solid line) together with the free spectrum (dashed line), (b) Profile in the space-line region with varying the chemical 
potential as n/K =0.59, 0.58, 0.575, 0.572, 0.571, 0.57006 (= /ic) (solid lines), 0.5682, 0.5652, 0.56, 0.55 (dashed lines). 



1 Tfl 

---.^/(.,0)=0. (17) 

In the chiral limit the energies of the stopped sigma and pion are w = 2M and 0, respectively, which fact is known as 
Nambu relation [19] in the symmetry-broken phase. 

At the CEP, on the other hand, the constituent quark mass M remains massive due to the non-zero current quark 
mass, and therefore we expect the sigma meson has a finite energy gap as a remnant of the Nambu relation. From the 
argument using the Landau free energy, M is expected to scale with m}/^ at the CEP [12,23]. With neglecting the 
imaginary part of /, we numerically confirm this behavior, in Fig. 2, for the sigma mass determined by Eq. (16) 
as well as the constituent quark mass M. The pion mass will be proportional to {m/My/^ ~ to^/^ from Eq. (17), 
and it is numerically found as ~ mP-'^^. It should be stressed here that the massive behavior of the sigma is quite 
normal from the viewpoint of the Nambu relation, and does not introduce any problem in the description of the phase 
transition. Actually it is found that the massive sigma and the vanishing curvature take place at the same time at 
the CEP in the NJL model [29], as well as in the linear sigma model [30], in the mean-field approximation. 



B. Spectrum of the scalar mode and the phase transition 

As we explained in Introduction, we expect that the instability occurs in the almost static and long wave length 
region in the scalar density fluctuation (as well as in the vector and entropy fluctuations through the coupling). In 
order to investigate the spectrum of the relevant modes in the transition at the CEP, we must keep the momentum 
q finite. 

The spectral function of the scalar mode, p{co,q) = 2ImxTOm('^) is explicitly written as 

1, , -(^)^(4Af^-g^)Im/(^,q) 

-Pico, q) - Imxmmiu;, q) - [_L „ j(q) + (4m2 _ g2)R,j(^, q)p + [(4^2 „ f)lmlic., q)]^ 

up to delta functions coming from possible non-decaying stable modes. We noticed previously that Ilmmi^,<l) is the 
scalar response function of the free quark gas of mass M, whose imaginary part 

^p^°\oj,(i) = Imn„„(a;,q) = -{AM^ - q^)lml{u,ci) (19) 

is the scalar channel spectrum of the uncorrelated quarks. These spectral functions have the support in the time-like 
region with = — q^ > 4M^ and in the space-like region with q"^ = w'^ — < 0. 

In the numerical calculation presented here, we fix the current quark mass to m/A = 0.01 for definiteness. The 
corresponding CEP locates at {T/A,p/A) = (0.14985,0.57006). The spectrum of XmrnC^^q) with |q|/A = 0.1 is 
shown in Fig. 3 (a), in which we see two sharp peaks at very low frequency and around 2M (M/A = 0.331). The free 
spectrum is also displayed for comparison. It is understood that the scalar channel attraction, which is encoded as 
a bubble sum in the denominator in Eq. (7) or in Eq. (18), changes the structure of the spectrum dramatically. The 
resulting scalar mode develops two distinct branches. The higher frequency peak corresponds to the sigma meson while 
the lower one appearing in the space-like region {w < \q\) is the scalar density fluctuating mode which we expect is 
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FIG. 4. (a) The pole position of Xmm(<^,q) as a function of |q|/A for the CEP a lower density {ii/h. = 0.56, x), and 
a higher density (/it/A = 0.58, +), together with the straight lines (oc |q|^ and |q|^). (b) Spectral function of the single pole 
form Eq. (21) (dashed line) in the space-like region at the CEP is compared with the numerical result (18) (solid line) with 
|q|/A = 0.1. 



responsible for the second order phase transition. Actually the latter peak diverges as q ^ 0. We study the structure 
of the peak in the space like region near the CEP in more detail in Fig. 3 (b). The enhancement of the spectrum is 
clearly demonstrated as the CEP is approached from the higher densities, /x/A = 0.59,0.58,0.575,0.572,0.571 (solid 
lines), and also from the lower /x/A = 0.55, 0.56, 0.5652, 0.5682 (dashed lines) with /Xc/A = 0.57006. 

We searched the pole in the unphysical lower half plane of lo [25-27], responsible for this peak. Indeed the pole is 
found on the negative imaginary axis as shown in Fig. 4 (a). We observed that the pole is moving toward the origin as 
decreasing the momentum q. Even in the non-critical case there is a pole on the negative imaginary axis and it moves 
toward the origin as q — * 0, too. The difference between the critical and non-critical cases exists in the exponent of 
the q-dependence; the pole position goes to zero as |q|^ in the critical case while on the other hand it behaves like 
|q|-^ for non-critical case. 

This result can be understood if one notes that in the static limit of the "sound" mode, a; ^ with u = a;/|q| < 1 
fixed, the function I has the form {v = k/E) [31] 



I{u) = 7(0, |q| -^0)-£^£ dk{n'_ + n'^)^ (in 



V + u 



V — u 



i7r6l(l - 



(20) 



whose imaginary part is proportional to u = w/|q| for small u. This imaginary part physically expresses the damping 
of the scalar "sound" mode in the quark gas. Then the response function with the space-like kinematics near the 
CEP may be approximated with 



Xmm(w,q) 



Renmm(0,q) 



1 



A(q) 



-i25(Imn„„(a;,q) -h (1 - 25(Ren„„(0,q)) _i^+;^„^(q) _it^ + t^;^(q) 



(21) 



where Xmm(q) is the scalar susceptibility, and Wc(q) = Xmm(fl)'^(fl)- Here we expect A(q) cx |q|, and numerically 
confirmed it. The susceptibility behaves as X^m(*l ^ 0) oc |q|^ at the CEP, but stays constant in non-critical region. 
This difference of the q-dcpendence of Xmm(q) alters the exponent of the q-dependence of the pole position, as 
described above and shown in Fig. 4 (a). 

We confirmed in Fig. 4 (b) that the w-dependence of the response function at the CEP is well reproduced in the 
single pole form (21) with the values (wcCq), Xmm(*l)) = (0.90925 • 10^^,0.21300 • 10~^) obtained in our numerical 
pole search. The deviation seen around uj/ K ^ 0.1 is due to the lack of the kincmatical cutoff to confine the spectrum 
inside the space-like region when we use the single pole form (21). We also checked that the response functions 
for the non-critical case near the CEP can be also described with the single pole form fairly well. Incidentally it is 
immediately seen that the form of (21) has the different limits around the origin: if one take the q — > limit first, 
XTOm(^j 0) 0. On the other hand, Xmm(0, q) = Xmm(q) ^ oo as q — > 0. 

The typical slow frequency of the system near the CEP is Wc(q) = Xmm('3)A('3)- This means that we reconfirmed the 
conventional theory of dynamical critical phenomena [32] in the form of Eq. (21): the divergence of the susceptibility 
is directly related to the critical slowing down of the response of the system to the external disturbance. The exponent 
of the dynamical scaling in our model will be equal to 3. 
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IV. SUMMARY 



In this paper we have explicitly shown that there are two branches in the scalar mode in the NJL model at finite 
(T, /i). The sigma meson mode stays massive at the CEP, following the approximate relation, Ad„ ^ 2M. The other 
mode is the scalar density fluctuation which becomes soft and relevant near the CEP. This result is perfectly consistent 
with the description of the second order phase transition using the Landau free energy. The first observation of the 
massive sigma at the CEP was done in Ref. [10] although it led the authors to a somewhat confused discussion. 

All the susceptibilities Xabici,b = T,ii,m,) arc found divergent at the CEP due to the flatness of the potential 
curvature, K ^ 0. The importance of the density fluctuation at the CEP is shown by considering the number- 
density- fixed susceptibility, Xmm{T, Pq). It no longer diverges at the CEP, but still does at the chiral phase transition. 
This situation at the CEP is quite similar to the liquid gas phase transition, where the density difference between 
two phases is idcntifled as the relevant order parameter. 

From our analysis we conclude that the CEP line in Fig. 1 denotes the second order phase transition characterized 
by the instability of the density fluctuations. Two lines of the CEP and the line of the second order chiral transition 
meet at the TCP (therefore tricritical), where we expect that the sigma meson and the density fluctuation mode 
become unstable at the same time. Although some readers could expect our conclusion from the beginning, there 
are confusions on this point in the literature [6,10,13] and one of our contributions is to clarify the importance of the 
density fluctuation explicitly using a concrete model for the CEP of a quark system. 

As for the experimental signatures of the CEP one should investigate first the implications of the anomalous 
fluctuations [7,8,12,14,33-35] of the baryon number, entropy, and scalar densities. One should be careful that the 
modification of the density fluctuations with the space-like momentum does not affect directly the spectrum of the 
particle production modes, like a tttt, II, 77, because of the kinematical mismatching. 
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